The generalized equation of motion for the orbital dynamics in the presence of current 
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The orbital dynamics induced by the charge current is studied theoretically. The equation of 
motion for the isospin vector T^ in the SU(N) case is derived in the presence of the current, and is 
applied to the cases of Cg (N=2) and t2g (N=3) orbitals. In spite of the anisotropic transfer integrals 
between orbitals, the dynamics is found to be isotropic for e^-orbitals similarly to the spin case , 
while it is anisotropic for t29-orbitals. The implication of this result to the current driven orbital 
domain wall motion is discussed. 



I. INTRODUCTION 

The dynamics of the spins in the presence 
of the current is an issue of recent intensive 
interesto^i^i^'^i^^i'^i^Qi^°i"i^^i". Especially it has been 
predicted theoreticallji^iiiSiiSiiiiiSii^ and experimentally 
confirmed that the magnetic domain wall (DW) motion 
is driven by the spin polarized current in the metallic 
ferromagnetsS and also magnetic semiconductors^^:. The 
basic mechanism of this current driven dynamics is the 
spin torque due to the current, which is shown to be 
related to the Berry phased. In addition to the spin, 
there is another internal degree of freedom, i.e., orbital 
in the strongly correlated electronic systems. Therefore 
it is natural and interesting to ask what is the current 
driven dynamics of the orbital, which we shall address 
in this paper. In the transition metal oxides, there 
often occurs the orbital ordering concomitant with the 
spin ordering. Especially it is related to the colossal 
magnetoresistance (CMR)M in manganese oxides. One 
can control this ordering by magnetic/electric field, 
and/or the light radiation. Even the current driven 
spin/orbital order melting has been observed. On the 
other hand, an orbital liquid phase has been proposed 
for the ferromagnetic metallic state of manganite&--. 
The effect of the current on this orbital liquid is also an 
interesting issue. 

There are a few essential differences between spin and 
orbital. For the doubly degenerate eg orbitals, the SU(2) 
pseudospin can be defined analogous to the spin. How- 
ever, the rotational symmetry is usually broken in this 
pseudospin space since the electronic transfer integral de- 
pends on the pair of the orbitals before and after the 
hopping, and also the spacial anisotropy affects as the 
pseudo-magnetic field. For t2g system, there are three de- 
generate orbitals, and hence we should define the SU(3) 
Gell-Mann matrices to represent its orbital state. There 
is also an anisotropy in this 8 dimensional order param- 
eter due to the same reason as mentioned above. 

In this paper, we derive the generic equation of motion 
of the SU(N) internal degrees of freedom in the presence 
of the orbital current. N=2 corresponds to the eg orbitals 
and N=3 to the t2g orbitals. Especially, the anisotropy 



of the DW dynamics is addressed. Surprisingly, there 
is no anisotropy for the eg case while there is for the t2g 
case. Based on this equation of motion, we study the DW 
dynamics between different orbital orderings. The effect 
of the current on the orbital liquid is also mentioned. 

II. CP"^^ FORMALISM 

In this section, we are only interested in the orbital 
dynamics. To be general, we consider a system with N 
electronic orbital degeneracy. The model that we inves- 
tigate is very similar to the lattice CP'^"^ sigma model 
with anisotropic coupling between the nearest-neighbor 
spinors. In contrary to the prediction in the band theory, 
the system is an insulator when it is 1/7V filled, namely 
one electron per unit cell. Because of the strong on- 
site repulsive interaction, the double occupancy is not 
allowed. Therefore, it costs high energy for electrons to 
move, and the spin degrees of freedom is quenched to 
form some spin ordering. Due to the complicate spin, or- 
bital, and charge interplay scheme, it is convenient to use 
the slave-fermion method in which we express the elec- 
tron as da-fi = h\z^^}z^^} where h\, z}^} , z^'^ are baptized 
as holon, spinon, and pseudo spinon (for the orbital) re- 
spectively, and the index i denotes the position^^ . If holes 
are introduced into the system, their mobility leans to 
frustrate the spin ordering and thus leads to a new phase 
which might possess finite conductivity. Therefore, we 
consider the following effective Lagrangian 

L = - hihi)ZaiZai -I- ^ {t"fhihjZaiZpj -I- c.c.)(l) 

i <ij> 

where z is short for the spinor z The iff in Eq.IH) 
is the transfer integral which is in general anisotropic 
because of the symmetry of the orbitals. Therefore, the 
system does not have SU(N) symmetry in general. To 
introduce the current, we consider the following mean 
field: 

< h,hj >= xe'^'J (2) 

where x denotes the doping concentration, and 9ij is the 
bond current with the relation 0ij = —Oji. Then, the 
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Lagrangian can be written as 

L = ih^{\ - x)Za^Za^ + ^ {t"f Xe'^'^ Z^^Zp^ + C.c).(3) 
i <ij> 

Note that the constraint J^a \^ai\'^ = 1 is imposed, and 
the lowest energy state within this constraint is reaHzed 
in the ground state. The most common state is the or- 
bital ordered state, which is described as the Bose con- 
densation of z, i.e., < Zai >=/= 0. In the present language, 
it corresponds to the gauge symmetry breaking. On the 
other hand, when the quantum and/or thermal fluctu- 
ation is enhanced by frustration etc., the orbital could 
remain disordered, i.e., the orbital liquid stated. Then 
the Lagrangian Eq. (|3Jl describes the liquid state without 
the gauge symmetry breaking. In this case, the gauge 
transformation 

^ai ^ 6 ^ ^ai (4) 

is allowed. Given 9ij — fij ■ j, where Vij — Vi — fj, 
the local gauge transformation in Eq. Q with ipi = ri ■ j 
corresponding to the simple shifts of the momentum from 
k to k + j. Therefore, the presence of current does not 
affect the state significantly since the effect is canceled 
by the gauge transformation. On the other hand, if Zai 
represent an orbital ordering, the current couples to the 
first order derivative of z in the continuum limit. Define 
a^t — izdfj^z, the second term in Eq.Q can be written 
as —j^a^ in the continuum limit. Namely, the current 
couples to the Berry's phase connection induced by the 
electron hopping. Therefore, we expect some non-trivial 
effect similar to the spin case. 

To derive the equation of motion for the orbital mo- 
ments, we use the SU(N) formalism. Introducing the 
SU(N) structure factors 

[A^,A^] = z/abcA^ (5) 
{A^,A^} = d^scA^-f5A<5^^ (6) 

where A^^ are the general SU(N) Gell-Mann matrices, [ ] 
are the commutators, and { } are the anti-commutators. 
The fABC is a totally-anti-symmetric tensor, and (Iabc 
is a totally-symmetric tensor. Let us express t"^^ in this 
basis 

U,^i^^l + t^)^ (7) 

We will only consider the nearest-neighbor hopping: 
tf, — i^,„-^. In the rest of the paper, , - will be written 

^tj^ ^ ^ i,i±k 

as t^, so does 9ij. Define the CP^^^ superspin vector as 

T^{i) = z^Xf^zp, (8) 

The equation of motion of T^{i) given by Eq.© can be 
obtained as^^. 

T''{^) = 7T^[-2cos0fe/wfT^(O 
(1 — xjn 

~ 2smektlAkT^{i) - sm9kt^dABcAkT^ii) 
+Wif sin0Uc(«)] (9) 



where the dummy k is summed over x, y, and z direction, 
a is the lattice constant, and the orbital current jc(*) is 
given by 

Jcii) = i{AZaiXapZ(3i - ZaiX'^pAzfji), (10) 

which is the second order in 9. Up to the first order in 9, 
the first term in the right hand side of the Eq.® is zero 
provided that 

t^+t^ + tf = (11) 

which is true for most of the systems that we are inter- 
ested in. Consequently, the dominant terms in Eq. (|5J| will 
be the second and the third ones, which can be simplified 
as 

T^{i) = --r^[2sm9,tlA,T^i^) 
(1 — x)n 

+ sm9kt^dABC^kT^{i)], 

(12) 

which is one of the main results in this paper. Using 
Eq. (|12() , we will discuss the orbital DW motion in the eg 
and the t2g systems. 

Here some remarks are in order on the mean field ap- 
proximation Eq.(|2Jl for the Lagrangian Eq.(^). First, it is 
noted that the generalized Landau-Lifshitz equation ob- 
tained by Bazaliy et. a& can be reproduced in the present 
mean field treatment when applied to the spin prob- 
lem. As is known, however, there are two mechanisms of 
current-induced domain wall motion in ferromagnetsii. 
One is the transfer of the spin torque and the other is 
the momentum transfer. The latter is due to the back- 
ward scattering of the electrons by the domain wall. Our 
present mean field treatment and that in Bazaliy's papes^ 
take the former spin torque effect correctly, while the lat- 
ter momentum transfer effect is dropped, since the scat- 
tering of electrons is not taken into account. However, 
the latter effect is usually small because the width of the 
domain wall is thicker than the lattice constant, and we 
can safely neglect it. 



III. N=2, eg SYSTEM 

First, we consider the application on the (La,Sr)MnO 
(113 or 214) system-'*. Without losing the generality, we 
show the result on the 113 system. The other one can be 
obtained in a similar way. 

In LaMnOs, the valence of Mn ion is Mn'^^ with the 
electronic configuration (t2g)^{eg)^ . By doping with Sr, 
one hole is introduced to Mn^+ and make it to be Mn^+ . 
The transfer integral between Mn ions depends on the Mn 
3d and O 2p orbitals. After integrating over the oxygen p 
orbitals, the effective hopping between the Mn d orbitals 
can be obtained. If we denote the up state as d^z^-r^ 
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and the down state as d^^-y^ , tij have the following form 



to(-l-— a - -a 




1 %/3 
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,1 , V3 , 1 
^2^ + -" ~t 



(13) 



where cr* are the Pauli matrices. For N=2, the pseudospin 
moment has the 0(3) symmetry given by 
The consequent equations of motion is given by 



/I 

za z. 



d X 



1 



(14) 



where 6 is {6^, dy,dz), which can be related to the orbital 
current as Vo = —j^^toO. Taking the continmmr limit, 
Ea. ((n|l becomes 



Moreover, due to the special properties of SU(2) alge- 
bra, the equation of motion is isotropic regardless how 
anisotropic the transfer integral is. Therefore, the mo- 
tion of orbital domain wall is undistorted as in the spin 
case. 




- Vo ■ A)T\r) = 



(15) 



which suggests the solution to be the form T'(r + Vot)- 
The result is similar to spin case. While the spin domain 
wall moves opposite to the spin current^^, in our case, 
the orbital domain wall also moves opposite to the orbital 
current. 

We can estimate the order of magnitude of the critical 
current to drive the orbital DW. The lattice constant a 
is about 3A. The transfer integral constant tg is around 
2eV in the LSMO system estimated from the photoemis- 
sion measurement. If we set v around Im/s, the critical 
current can be estimated as e/a'^ ~ 6 x 10^ A/m^, which 
is roughly the same as the order of magnitude of that to 
drive the spin domain wall. 

It should be noted that the current only couples to 
the first order derivative of z. The double exchange term 
which is given by the second order derivative is not shown 
in the equation of motion. However, the double exchange 
term plays a role to stabilize the DW configuration before 
the current is switched on. 

There are two orbitals degenerate to dx^-y^, which are 
rfy2„^2 and (iz2_j.2. Similarly, d3y2_^2 and d^^i^^i are 
degenerate to d^z^_^2. In the Manganite system, d^z2_j.2 
and dx2_y2 may have different energy due to slight struc- 
tural distortion in the unit cell. Therefore, in most cases, 
domain walls are the type of those which separate two de- 
generate domains. For example, let's consider the orbital 



domain wall separating iy^ — and 3x^ 



In FigCl 



the domain wall sits at ?/ = 0, and suppose the current 



is along the positive x direction. 3y 



and 3x — 



orbitals are described by the spinors (—1/2, — \/3/2) and 
(— 1/2, -\/3/2) respectively. The configuration given in 
Figgis described by the spinor field (cos6'(xi), s\n9{xi)) 

1), where w is the 



where 9{xi) = 27r/3(2cot 
width of the domain wall. 



-Xi j w 



FIG. 1: (Color online) An example of the domain structure 
in the eg system. The domain boundary lies on the y-axis. 



IV. N=3, t2a SYSTEM 

Let us consider the tig systems, for example, in the 
Vanadate or Titanate systems. The tig orbitals contain 
three orbits d^y, dyz-, and d^x- The hopping integral t^ 
between the Ti'^+ sites or the V'^+ ones is given as 



tx, — 




= to(2i + ^A3-^A8) 
^3 2 2a/3 ' 



^M^1 + ^A«) 



^3" 2 ^ 2Vi' 



= t,C-l-^X'--^\^) 



(16) 



where A'^ are SU(3) Gell-Mann matrices with the nor- 
malization condition Tr (A'^A^) = 25^^ . The super-spin 
T"* given by zX^z is normalized because zz = 1 and 

The equation of motion in this case will be anisotropic 
because dABC is non-trivial in the SU(3) case. It is in- 
spiring to work on an example to see how it goes. Let's 
consider the dxy — dyz orbital DW shown in the Fig|21 
Because of the orbital symmetry, such DW can only be 
stabilized along the y— direction. Similarly, the dxy — d^x 
DW can only be stabilized along the x— direction and so 
on. Therefore, the effect of current is anisotropic. In the 
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absence of current, the domain wall is stabilized to be 



sin tan 



Za{ri) — I costan y^^'" 




(17) 



It can be easily seen that it takes no effect if the current is 
applied along a; or z direction. The superspin component 
is given by 

/ sech(y^/w) \ 


tanh(yj/w) 


T^(f) = (18) 






1 
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When the current is applied along the y— direction, 
Eq.jnj for e; 
are given by 



Eq.(|T21) for each are decoupled. For A = 1,2, 3, they 



(7^ + 



2xtQ 



d 



dt {l-x)h "dy 



)Ti'2'3(f) = 



For A = 4. .8, they are given by 



xto 



d 



(— 



)T^-^{r) = 



(19) 



(20) 



At a glance, we obtain 2 characteristic drift velocity of 
the domain wall. It is not the case, because r^(r) is 
zero and T^{r) is constant. Only Ea. H19|l determines 
the motion of the DW . Furthermore, the 8— dimensional 
super-spin space reduces to be 2— dimensional as sum- 
marized in FigI3 moment grows in the domain wall 
while moment distinguishes two domains. In the 
presence of current along direction, the wall veloc- 
ity is \v\ = {i-^x)h \^y\- Other type of domain struc- 
tures can be analyzed in a similar way. As a result, even 
though the order parameter in the t2g systems forms an 
8-dimensional super-spin space, we can always reduce it 
to be 2-dimensional because of the anisotropic nature of 
the system. Furthermore, the DW moves without any 
distortion just like the isotropic case. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we formulated the orbital dynamics when 
the spin degree of freedom is quenched. We used SU(N) 
super-spin to describe the orbital states and obtained 
the general equation of motion for it. We also showed 
some examples for the SU(2) and SU(3) cases correspond- 
ing eg and t2g systems, respectively. In the SU(2) case, 
the orbital dynamics is very similar to the spin case: 
undistorted and isotropic. In SU(3) case, the DW struc- 
ture is anisotropic because of the orbital symmetry. In 




FIG. 2: (Color online) An example of the domain structure 
in the t2g system. Only y— component of current will move 
the domain wall. 




FIG. 3: (Color online) dxy — dyz domain wall in the superspin 
space. The superspin rotates like a XY spin in the superspin 
space. 



addition, the effective super-spin space is 2-dimensional, 
and the domain wall motion is also undistorted. 

Even though the analogy to the spin case can be made, 
one must be careful about some crucial differences be- 
tween the spin and orbital degrees of freedom. We have 
estimated the critical current to drive the domain wall 
assuming the uniform current flow in the metallic sys- 
tem, but most of the orbital ordered state is insulating. 
This is the most severe restriction when the present the- 
ory is applied to the real systems. An example of the 
metallic orbital ordered state is A-type antiferromagneitc 
state with — orbital ordering in NdSrMnOii. How- 
ever it is insulating along the c-direction, and there is no 
degeneracy of the orbitals once the lattice distortion is 
stabilized. The ferromagnetic metallic state in LSMO is 
orbital disordered. According to the quantum orbital liq- 
uid picture -^, there is no remarkable current effect on the 
orbitals as explained in the Introduction. On the other 
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hand, when the classical fluctuation of the orbital plays 
the dominant role for the orbital disordering, the short 
range orbital order can be regarded as the distribution of 
the domain walls, which shows the translational motion 
due to the current as discussed in this paper. Note how- 
ever that the current does not induce any anisotropy in 
the orbital pseudospin space in the Cg case. 

Orbital degrees of freedom is not preserved in the vac- 
uum or the usual metals, where the correlation effect is 
not important, in sharp contrast to the spin. There- 
fore the orbital quantum number can not be transmitted 
along the long distance and the pseudospin valve phe- 
nomenon is unlikely in the orbital case. 
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APPENDIX A: ORBITALS 



There are two class of Cg orbitals, i.e., one is ix^ — r^, 
3y2 _ j,2^ _ j,2^ g^j^j ^Yie other is — j/^, ~ z^ , 

and z^ — . They can be described by a two-component 



spinor. If we use (1,0) to describe 2>z 



and (0, 1) for 



y , the other degenerate orbitals are given by 

3y2-r2 : (-1/2,-^3/2) 

Zx^-r^ ■ (-l/2,\/3/2) 

v'-z' : (-V3/2,-l/2) 

: (V3/2,-l/2) (Al) 



The Gell-Mann matrices used in the text are given ex- 
plicitly in the following: 



Ai = 



A3 = 




• (Bl) 
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